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SECTION  I 


INTRODUCTION 

Dramatic  increases  in  the  size  and  flexibility  of  new  aircraft  have 
induced  a  greater  reliance  upon  power  spectral  methods  in  gust  design. 
Corresponding  structural  representations  of  the  aircraft  are  typically  of 
the  lumped  mass  variety,  and  three-dimensional  lifting  surface  aero¬ 
dynamics  is  routinely  employed.  However,  in  certain  important  respects, 
the  current  atmospheric  turbulence  description  has  not  kept  pace  with 
these  developments. 

Since  gust  velocity  varies  in  magnitude  and  direction  along  the  path  of 
an  airplane,  it  is  obvious  that  a  similar  variation  must  exist  along  any  line 
perpendicular  to  the  path.  It  has  long  been  recognized  that  the  latter  gust 
variation  may  be  sufficient  over  the  transverse  dimensions  of  a  large 
aircraft  to  justify  its  inclusion  in  the  analytical  turbulence  description. 
However,  because  of  the  extremely  cumbersome  mathematics  required  to 
generate  the  necessary  turbulence  description,  the  most  common  procedure 
has  consisted  of  retaining  only  the  vertical  gust  component  and  ignoring  gust 
variations  which  depend  upon  vertical  position.  This  results  in  a  single¬ 
component,  two-dimensional  turbulence  model.  The  earlier  and  more 
common  one-dimensional  model  completely  ignores  transverse  variations 
in  gust  and  treats  vertical  and  lateral  gust  responses  as  separate  calcula¬ 
tions.  ! 

An  exact,  but  simple  three-dimensional  gust  model  consistent  with 
current  structural  and  aerodynamic  computational  procedures  is  presented 
in  References  1  through  4.  The  original  development  is  contained  in 
References  1  and  2,  and  consists  of  the  following  essential  steps:  (1)  a 
unique  rotational  transformation  is  applied  to  reduce  from  nine  to  four, 
the  number  of  independent  components  required  to  specify  the  gust  velo¬ 
city  cross  spectrum  tensor.  (2)  Closed  form  solutions  of  the  cross 
spectrum  tensor  are  derived  for  both  the  Dryden  and  von  Karman  spectral 
models.  (3)  The  computational  redundancies  which  arise  when  the  analytical 


1 


representation  of  the  aircraft  is  introduced  are  eliminated  by  invoking  the 
appropriate  symmetry  properties  of  the  cross  spectrum  tensor. 

This  procedure  results  in  a  concise,  unified  multiple  input  formulation  that 
describes  the  dynamic  response  of  a  bilaterally  symmetric,  flexible  aircraft 
traversing  a  random,  isotropic  field  of  atmospheric  turbulence.  The  new 
method  incorporates  in  its  present  form  the  following  significant  features; 

(1)  The  formulation  is  fully  three-dimensional,  so  that  even  if  the  aircraft 
is  assumed  to  be  unresponsive  to  one  or  two  of  the  three  components  of 
turbulence,  the  result  still  differs  from  the  1-dor  2-d  case.  (2)  The  com¬ 
putation  is  efficient;  i.e.,  all  three  components  of  turbulence  may  be  included 
in  a  single  calculation,  and  the  cost  per  case  of  performing  a  3-d  gust  response 
analysis  of  a  flexible  airplane  ie  essentially  the  same  as  the  corresponding  1-d 
calculation.  (3)  The  correlation  properties  of  the  3-d  turbulence  field  are 
conveniently  described  in  terms  of  only  two  functions  of  a  single  variable. 

These  are  designated  as  the  planar  and  nonplanar  gust  coherences  and  are 
furnished  in  tabular  form.  (4)  Since  the  formulation  separates  the  amplitude 
properties  of  the  turbulence  field  from  its  correlation  properties,  the  same 
analytical  or  empirical  gust  spectra  may  be  employed  as  in  the  conventional 
1-d  analysis,  (5)  If  the  scale  of  turbulence  is  much  larger  than  the  transverse 
dimensions  of  the  aircraft,  then  the  formulation  becomes  independent  of  the 
scale  except  in  specifying  the  analytical  gust  spectra.  (6)  If  the  transverse 
dimensions  of  the  aircraft  are  ignored,  then  the  method  reduces  to  the  con¬ 
ventional  1-d  analysis.  (7)  If  the  vertical  dimensions  of  the  aircraft  are 
ignored  and  the  calculation  is  confined  to  vertical  gust  only,  then  the  method 
reduces  to  the  current  2-d  analysis.  (8)  The  method  is  compatible  with 
current  turbulence  design  criteria  procedures,  so  that  only  a  possible  adjust¬ 
ment  of  existing  parameter  values  would  be  required.  (9)  Anisotropic  turbu¬ 
lence  can  be  introduced  by  a  simple  procedure  which  preserves  isotropy  at 
small  wavelengths. 

The  theory  encompasses  two  types  of  response  analysis:  power  spectral  and 
cross  spectral.  The  power  spectral  approach  may  be  used  to  furnish  data 


for  gust  design,  including  fatigue,  design  loads,  stability  and  control. 

Cross  spectral  methods  are  primarily  employed  to  relate  structural 
responses  to  gust  inputs  as  measured  by  a  probe  attached  to  the  aircraft. 

By  comparing  cross  spectral  data  with  actual  flight  measurements,  a 
direct  verification  of  the  analytical  model  may  be  achieved.  Power 
spectral  data  may  also  be  used  for  this  purpose.  The  work  described 
herein  covers  both  methods,  and  the  objective  is  to  demonstrate  and 
evaluate  the  effect  of  three-dimensional  turbulence  models  in  comparison 
to  one -dimensional  models  by  means  of  dynamic  response  test  data  from 
the  C-5A  aircraft. 

This  effort  represents  part  of  a  long  range  plan  which  is  programmed  to 
achieve  four  interrelated  objectives:  (1)  implementation  of  three-dimen¬ 
sional  gust  response  methods,  (2)  validation  of  these  methods  by  comparison 
with  dynamic  response  test  results,  (3)  incorporation  of  these  methods  into 
current  design  procedures,  (4)  development  of  new  methods  and  applications, 
including  treatment  of  aircraft  response  to  turbulence  conditions  which 
depart  from  current  ideal  assumptions.  This  report  is  intended  to  achieve 
objective  (2),  to  demonstrate  and  evaluate  the  effect  of  three-dimensional 
turbulence  by  means  of  available  aircraft  dynamic  response  data. 


SECTION  n 


GUST  VELOCITY  CROSS  SPECTRUM 

The  turbulence  description  required  for  a  3-d  gust  response  analy  sis 
of  an  aircraft  is  provided  by  the  gust  velocity  cross  spectrum  tensor  as 
measured  by  traversing  a  statistically  stationary,  random,  isotropic  field 
of  atmospheric  turbulence.  The  cross  spectrum  tensor  is  used  to  compute 
the  gust  velocity  cross  spectrum  between  gusts  impinging  upon  any  two 
points  of  the  aircraft.  The  analytical  representation  of  the  gust  velocity 
cross  spectrum  tensor  for  the  von  Karman  turbulence  model  was  originally 
derived  in  Reference  2  (See  also  References  3  and  4).  A  derivation  of  the 
cross  spectrum  tensor  for  the  Dryden  model  originally  appeared  in  Refer¬ 
ence  1,  but  it  is  considered  to  be  less  representative  of  atmospheric 
turbulence  under  most  flight  conditions. 

2.  1  CROSS  SPECTRUM  TENSOR 

A  unique  rotational  transformation  was  also  employed  in  References 
1  through  4  to  reduce  the  required  number  of  distinct  cross  spectrum  com¬ 
ponents  to  four.  Moreover,  the  transverse  separation  in  the  original  frame 
becomes  the  lateral  separation  in  the  rotated  frame,  so  that  explicit  depen¬ 
dence  upon  the  vertical  separation  variable  is  eliminated.  It  was  also  shown 
that  under  Taylor's  "frozen  turbulence"  hypothesis,  the  dependence  of  the 
cross  spectrum  tensor  upon  the  longitudinal  separation  reduces  to  the  con¬ 
ventional  phase  factor  for  gust  penetration.  As  such,  the  latter  quantity  may 
be  transferred  to  the  frequency  response  function  of  the  aircraft.  The  four 
distinct  cross  spectrum  components  are  thus  reduced  to  functions  of  the 

/~ 2  2 

reduced  frequency  fl  =  w/V  and  the  transverse  separation  p  ==  /y  +  z  . 

The  cross  spectrum  between  gust  velocity  components  measured  along 
the  respective  unit  directional  vectors  m  =  anc*  i!  =  ^ni,n2,n3^ 

which  are  separated  by  the  vector  r  =  (x,y,z)  whose  transverse  component 
is  p  =  (0,y,z)  ,  then  becomes 
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where 


4*  (p ,  f2)  =  *(e^e^,pe2rfl) 


(2.2) 


is  the  cross  spectrum  between  cartesian  gust  velocity  components  having  a 
purely  lateral  separation.  The  cartesian  unit  vector  is  denoted  by  e..  . 
The  vector  configuration  is  illustrated  by  the  following  sketch. 


-1 

ORIGINAL  FRAME  ROTATED  FRAME 
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2.  2  CLOSED  FORM  SOLUTION 


For  evaluation  and  analysis,  it  is  convenient  to  separate  the  correla¬ 
tion  and  amplitude  properties  of  the  four  distinct  cross  spectrum  compon¬ 
ents  of  the  tensor  by  factoring  them  into  the  general  form 


where 


and 


*jk(p,a)  =  $j]c(p,n>/*j(a)*k(n) 
$jk(P,fl)  =  *jk(pfn)//*j(n)*k(n) 

* .  («)  =  *  .  .  (0,fl) 

3  33 


(2.2) 

(2.4) 

(2.5) 


Equation  (2.  5)  defines  the  power  spectrum  of  a  gust  velocity  (the  auto 
spectrum  of  a  velocity  component).  It  is  a  real  function  representing  the 
frequency  distribution  of  the  mean  square  amplitude  of  a  cartesian  velocity 
component. 

Equation  (2.4)  defines  the  gust  coherence.  It  represents  the  normalized 
form  of  the  corresponding  gust  velocity  cross  spectrum  of  equation  (2.3), 
and  is  obtained  by  effectively  dividing  out  the  frequency  distributions  of  the 
respective  RMS  gust  velocity  components.  What  remains  is  a  complex 
function  whose  modulus  and  phase  angle  correspond  to  the  frequency  distri¬ 
butions  of  the  correlation  coefficient  and  phase  difference,  respectively, 
between  the  gust  components. 

References  2  through  5  give  the  following  closed  form  solution  for  the 
von  Karman  power  spectra  and  gust  coherences  required  in  equation  (2.  3) 
to  compute  the  four  cross  spectrum  components  of  equation  (2.  1). 
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1 


Vn)  " 


2La^/7T(l  +  <2)5/  6 


(j  =  1) 


1  La^  ( 1  +  8k2/3)  /  tt  ( 1  +  k2)11/6  (j  =  2,3) 


(2.6) 


whe  re 


K  =  1. 339LQ,  R  =  w/V 


and 


1.339  =  T (l/3)/r (1/2) T  (5/6) 


Also,  $i;L(p,R)  =  0.9944(p5/6  K5/g  (y)  -  y11/6  K1/g(y)/2) 


-i$12(p,R)  =  0 . 9944y11/  6  K5/g  (y)  </21/2/l  +  8k2/3 


(2.7) 


*22(p,R)-*33(p,R)  =  0.9944y11/6  K1/g(y)(l  +  +  8k*/3) 


$33(p,ft)  =  0.9944  (y5/6  K5/g(y)  -  y11/6  K1/g(y)/(l  +  8k2/3)) 


where 


y  =  (p/1. 339L) 


and 


0.9944  -  21/6  T(5/6) 


r  and  K  denote  the  gamma  function  and  the  modified  Bessel  function 
of  the  second  kind,  respectively.  Notice  that  the  scale  of  turbulence  L  , 
which  is  a  measure  of  average  eddy  size,  emerges  as  a  parameter  in  the 
closed  form  solution.  Since  the  coherence  between  longitudinal  and  lateral 
gust  is  a  positive  imaginary  quantity,  it  is  multiplied  by  -i  to  yield  a 
positive  real  quantity.  Similarly,  the  coherence  between  lateral  components 
is  replaced  by  a  difference  for  computational  convenience. 
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Figure  1  contains  the  longitudinal  and  transverse  (lateral  and  vertical) 
power  spectra.  They  are  normalized  and  plotted  versus  nondimens ional 
reduced  frequency  flL  ,  by  taking  the  RMS  gust  velocity  component  and  the 
scale  of  turbulence  to  be  of  unit  magnitude. 

Figures  2  and  3  show  the  four  gust  coherences  plotted  versus  pft  , 
the  lateral  separation  in  radians.  Each  coherence  appears  as  a  family 
of  curves  with  p/L ,  the  lateral  separation  in  scale  lengths,  as  parameter. 
All  of  coherences  are  monotonic  with  respect  to  the  parameter  except  for 
the  family  of  curves  at  the  top  of  Figure  3.  There,  for  values  of  the  para¬ 
meter  greater  than  are  shown,  the  curves  fall  back  onto  the  horizontal  axis. 

A 

The  fact  that  is  imaginary  indicates  that  a  tt/2  phase  differ¬ 

ence  exists  between  the  coherent  portions  of  the  longitudinal  and  lateral  gust 
velocity  components.  Similarly,  drops  below  the  zero  axis  and 

exhibits  small  negative  values  at  all  subsequent  lateral  separations.  There¬ 
fore,  coherent  vertical  gust  components  have  a  phase  difference  of  tt  under 
these  conditions.  Reference  to  equation  (2.  1)  shows  that  introducing  a  longi¬ 
tudinal  separation  merely  imposes  a  further  difference  in  phase  due  to  gust 
penetration,  and  has  no  effect  upon  the  modulus  of  the  coherence.  Moreover, 
since  the  modulus  corresponds  to  a  correlation  coefficient,  its  value  cannot 
exceed  unity. 

Reference  4  shows  that  although  the  power  spectra  for  the  Dryden  and 
von  Karman  spectral  models  are  quite  dissimilar,  the  numerical  values  of 
their  gust  coherences  are  almost  identical.  This  suggests  that  the  gust 
coherences  of  equation  (2.7)  may  be  used  in  equation  (2.3)  with  empirical 
power  spectra  that  need  not  conform  to  the  von  Karman  model. 

2.  3  LARGE  SCALE  LIMIT 

If  the  scale  of  turbulence  is  large  compared  to  the  transverse  dimen¬ 
sions  of  the  aircraft,  then  the  gust  coherences  of  equations  (2.7)  lose  their 
dependence  upon  the  scale  and  reduce  to  the  following  single -variable 
functions  of  the  lateral  separation  in  radians: 
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Figure  Z.  Coherences  of  Longitudinal  and  Lateral  Gust 
Components  for  the  von  Karman  Case 
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1 


*n(p*l)  =  0.9944  ((pfl) 5/6  K5/6(pO)  -  (pH) 11/6  K1/fi(pfi)/2] 


1/6 


-ii12<pft)  =  0.9944  (/J/4)  (pfl) 11/6 


K5/6(pG) 


(2.8) 


11/6 


>22  (p^2) -^33  (p«)  =  0.9944  (3/8)  (pft)  '  K  (pft) 


$33(pn)  =  0.9944  (pfi)5/6  K5/6(pQ) 


where  $.R(p^)  =  lim  $  .k(p,Q)  =  lim  $.,(p,f2) 

-*  p/L^O  ^  L-*-00  ^ 


(2.9) 


Notice  that  the  same  result  may  be  derived  from  the  more  restric¬ 
tive  assumption  of  infinite  scale.  It  is  therefore  referred  to  as  the  large 
scale  limit. 

2.  4  PLANAR  AND  NONPLANAR  COHERENCES 


Examination  of  the  four  terms  in  the  cross  spectrum  formulation  given 
by  equation  (2.  1)  shows  that  if  longitudinal  gust  is  ignored,  the  first  two 
terms  may  be  omitted.  Furthermore,  when  the  resulting  gusts  are  measured 
normal  to  a  plane,  the  fourth  term  assumes  its  maximum  value,  and  the  third 
term  vanishes.  We  will  therefore  designate  the  corresponding  cross  spectra 
as  the  planar  and  nonplanar,  respectively,  denoted  by  $p  and  .  Accord¬ 
ing  to  equations  (2.8)  and  (2.9),  the  corresponding  planar  and  nonplanar 
coherences  in  the  large  scale  limit  may  be  written 


*p(pn)  =  4>33  (pQ)  =  0.9944  (pfl)  5/6  K5/6(pft) 

*N(pfl)  =  $22  (pf2)  -  $33(pfi)  =  0.9944  (3/8)  (pfi)11/6  K1/6  (pfl) 


(2.  10) 
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The  asymptotic  coherences  defined  by  equations  (Z.  10)  may  be  sub¬ 
stituted  back  into  equations  (2.7)  by  changing  the  variable  from  pH  to  y 
and  matching  up  similar  terms.  All  four  coherences  are  thereby  expressed 
in  terms  of  the  asymptotic  planar  and  nonplanar  coherences.  The  result  is 


$11(p,fi) 


-i$12 (P/ft) 


$22  (p,fl) -$33 (p,Q) 
$33 (p,Q) 


=  *p(y)  -  (4/3)*N(y) 

=  *  (y)  \k/21/2A  +  8<2/3 
=  *N(y) (1  +  <2 ) / (3/8  +  k2) 
=  *p(y)  -  (p)  /  (3/8  +  k2) 


(2.11) 


Notice  that  4>p(pft)  and  $>N  (pft)  correspond  to  the  p/L  =  0  curves  in 
Figure  3.  As  single -variable  functions,  they  may  be  conveniently  tabulated 
for  use  in  the  above -express ions.  Furthermore,  if  p/L  is  sufficiently 
small,  then  the  gust  coherences  may  be  obtained  directly  from  equations 
(2.8). 


2.  5  SMALL  AIRCRAFT  LIMIT 

If  the  large  scale  limit  described  in  Paragraph  2.3  is  satisfied,  and 
if  the  transverse  separation  in  radians  is  sufficiently  small,  then 


lim  $>..  (p,fl)  =  lim  (p,fi) 
P/L-*-0  3K  p+0  DK 

pft-*-0 


vo) 


'  X 

.  0 


j  =  k 
j  /  k 


(2.  12) 


Notice  that  the  same  result  may  be  derived  from  the  more  restrictive 
assumption  of  zero  transverse  separation.  It  is  therefore  referred  to  as 
the  small  aircraft  limit,  although  the  restriction  applies  only  to  the  trans¬ 
verse  dimensions  of  the  aircraft  and  not  to  its  length.  Since  ft  =  co/V  ,  the 
limit  pft  =  0  is  also  approached  when  the  radial  frequency  is  small  or  the 
airspeed  is  large. 

2.  6  0-D,  1-D  AND  2-D  APPROXIMATIONS 

The  zero-dimensional  approximation  for  the  gust  velocity  cross  spectrum 
corresponds  to  the  assumption  of  a  point  aircraft  and  is  valid  if  gust  varia¬ 
tions  are  negligible  within  the  spatial  boundaries  of  the  aircraft.  Setting 
r  =  0  in  equation  (2.  1)  and  applying  equation  (2.  3)  and  (2.  12)  yields 

3 

4>{m,n,0^ft)  =  l  $.(ft)m.n.  (2.13) 

j=l  J  J  J 

(0-D  TURBULENCE) 


The  conventional  1-d  turbulence  approximation  corresponds  to  the 
assumption  of  a  linear  aircraft  and  is  valid  if  gust  variations  are  negligible 
except  across  the  longitudinal  dimension  of  the  aircraft.  This  condition  is 
identical  to  the  small  aircraft  limit  discussed  in  Paragraph  2,5.  Setting 
r  =  xe^  in  equation  (2.  1)  and  proceeding  as  before,  we  obtain 


0 (m,n,xe^,fl) 


-iftx 

e 


3 

l 


j=l 


$  .  (ft) m .n  . 
3  3  3 


(2.  14) 


(1-D  TURBULENCE) 


This  result  differs  from  the  O-d  turbulence  case  only  in  that  the  phase 
factor  for  gust  penetration  is  now  included. 

The  conventional  2-d  turbulence  approximation  corresponds  to  a  planar 
aircraft  and  is  valid  if  gust  variations  across  the  vertical  dimension  of 
the  aircraft  are  small  enough  to  be  ignored.  Setting  r  =  xe^  +  ye2  in 
equation  (2.1)  yields 


$  (m ,  n , xe ^+ye2 , ft ) 


=  e-^^x ($11  (y ,ft)m^n^  +  $12(y,ft)  (m1n1  +  m2n2)  + 

^22(y,fi^m2n2  +  $33(y»fi)m3n3)  (2.15) 

(2-  D  TURBULENCE) 


A  comparison  of  equations  (2.  1)  and  (2.  13-2.  15)  shows  that  the  0-d  and 
1-d  cases  require  only  the  gust  power  spectra,  whereas  the  2-d  and  3-d  cases 
also  require  the  gust  coherences  according  to  equation  (2.  3). 

2.  7  OMISSION  OF  GUST  COMPONENTS 

The  aerodynamic  representation  of  the  aircraft  to  which  the  gust  velocity 
cross  spectrum  is  applied  may  exclude  response  to  a  given  cartesian  gust 
velocity  component.  In  that  case  the  corresponding  components  of  the  direc¬ 
tional  vectors  m  and  n  in  the  cross  spectrum  expression  are  simply  set  to 
zero.  This  may  result  in  the  elimination  of  various  terms  .in  equations  (2. 1) 
and  (2.13-2.15).  For  example,  in  Reference  5  the  2-d  approximation  is 
applied  to  the  Concord  aircraft,  but  longitudinal  and  lateral  gust  responses 
are  omitted.  In  that  case  equation  (2. 15)  reduces  to 


B  4  Q  y 

$  (m,n,xe^+ye2  ,ft)  =  e  4>(y,ft)m3n^ 
(2-D  VERTICAL  GUST) 


(2.  16) 
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The  response  to  lateral  gust  may  be  included  in  a  separate  1-d  calculation. 
In  that  case  equation  (2.  14)  reduces  to 


t(mfn,xeirn)  =  e_lfix  $2(ft)m2n2 

(1-D  LATERAL  GUST) 


(2.  17) 


2.8  ANISOTROPIC  TURBULENCE 

Both  the  scale  of  turbulence  and  the  RMS  value  of  a  gust  velocity  com¬ 
ponent  may  depend  upon  direction  if  the  turbulence  field  is  anisotropic.  In 
that  case  the  following  substitutions  may  be  made  in  the  von  Karman  power 
spectra  for  isotropic  turbulence  given  by  equations  (2.6). 


°W  -  °j  '  L  *  Lj 


(ANISOTROPIC  TURBULENCE) 


(2.  18) 


Constraints  upon  the  values  of  these  parameters  may  be  derived  from 
the  condition  that  anisotropic  turbulence  approaches  ;sotropy  in  the  high 
frequency  range  (Reference  6).  Therefore,  the  gust  spectra  for  anisotropic 
turbulence  must  assume  the  same  asymptotic  form  as  equations  (2.6)  in  the 
limit  of  large  ft,  thereby  yielding 


L  . 
J 


«  a3 
3 


(2. 19) 


Moreover,  if  horizontal  isotropy  is  present,  then  the  further  restriction  that 
=  L2  may  be  imposed.  These  constraints  may  be  invoked  when  applying 
an  analytical  fit  to  simultaneously  measured  gust  spectra.  They  are  also 
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F* 


useful  for  establishing  gust  criteria  for  analytical  gust  spectrum  formula¬ 
tions  and  are  reflected  in  current  military  design  criteria  specifications 

(References  7  and  8). 


Figures  2  and  3  show  the  gust  coherences  to  be  quite  insensitive  to  the 
assumed  scale  of  turbulence  for  the  small  values  of  p/L  that  typically  pre¬ 
vail.  It  is  then  permissible  to  assume  the  isotropic  turbulence  relation 
L  =  in  computing  the  gust  coherences  for  anisotropic  turbulence. 
Anisotropy  in  gust  amplitude  may  be  accounted  for  by  multiplying  each 
directional  vector  component  nr  and  by  the  corresponding  power  spectrum 


amplitude  /4>j  (ft)  and  in  compensation,  replacing  each  cross  spectrum 


$.  (p,ft)  by  its  coherence  4|'t(p,ft). 

3  k  jk 


i 

) 
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SECTION  HI 


AIRCRAFT  FREQUENCY  RESPONSE  FUNCTIONS 

The  aircraft  structural  and  aerodynamic  description  used  in  the  analy¬ 
sis  of  flight  in  turbulence  is  contained  in  the  conventional  frequency  response 
function  which  relates  a  given  structural  response  to  the  gust  impinging  upon 
an  aerodynamic  reference  panel.  This  function  is  independent  of  the  turbulence 
description  and  is  therefore  applicable  to  all  of  the  turbulence  cases  discussed 
in  Section  II.  The  necessary  frequency  response  functions  for  loads  and 
accelerations  are  derived  in  the  following  paragraphs  from  the  principle  of 
balanced  forces.  Matrix  notation  is  used  throughout  the  development. 

3.  1  FORCE  EQUILIBRIUM 

The  aerodynamic  surfaces  of  the  aircraft  are  assumed  to  be  represented 
by  reference  panels  with  respect  to  which  there  exists  a  normal  gust  velocity 
distribution  {w}  that  is  furnished  independently  by  the  turbulence  description. 
The  resulting  generalized  displacements  {£}  are  assigned  to  the  mo',n  °hapes 
derived  from  a  free-free  modal  analysis  of  the  aircraft  structure.  The  force 
balance  principle  for  the  aircraft  may  be  expressed  in  terms  of  either  the 
gene ralized  forces. 


(q/V)  [Gj  fe  inx>J{W}  + 

(oj2  [M]  -  tMjrw2(l  +  ig)J  +  q  [  A]  +  q[C]  [F]  [S])  (C)  =  (0}  (3.1) 


or  the  response  loads 


(q/V)  [G]  re‘ifix>J{w}  + 

(w2[M]  -  [M]fiJ2(l  +  ig)J  +  q  (A)  +  q  [C]  [F]  [S]  )  {£  }  =  {0}  (3.2) 
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The  first  term  of  each  equation  represents  the  applied  gust  forces,  whereas 
the  second  term  contains  the  forces  resulting  from  the  displacement  of  the 
structure.  These  are  obtained  by  multiplying  the  generalized  displacements 
by  the  respective  inertia,  elastic,  structural  damping,  aerodynamic  and  con¬ 
trol  surface  load  coefficients.  Notice  that  the  load  coefficients  which  furnish 
the  generalized  forces  are  indicated  by  an  overbar  in  equation  (3.  1)  to  distin¬ 
guish  them  from  those  which  yield  the  response  loads  in  equation  (3.  2). 

Noti'.e  also  that  a  diagonal  matrix  of  gust  penetration  factors  has  been  inserted 
to  shift  the  phase  of  {w}  from  the  origin  of  the  aircraft  coordinate  system  to 
the  aerodynamic  reference  panels. 

3. 2  ACCELERATIONS 

Equation  (3.  1)  may  be  solved  for  the  generalized  displacements  by  matrix 
inversion.  Thus, 


(U  =  (q/v)  [S'] -1  [G] 


(3.3) 


where  [B]  =  —  fMj  (oj2  -  [“w2(l  +  ig)J  )  -  q  [A]  -  q[C][F][S] 


The  matrix  [5]  ^  is  seen  to  contain  the  ratios  of  generalized  displacement 
to  generalized  gust  force. 

Let  [R]  represent  the  ratios  of  accelerometer  displacement  to  general¬ 
ized  displacement.  Then  with  the  aid  of  equation  (3.  3),  the  acceleration 
responses  may  be  written 


-to2  [R]  {£}  =  -U)2  (q/V)  [R]  B]_1[G]  (3.4) 


The  matrix  coefficient  of  (W}  in  the  above  equation  contains  the  influence 
coefficients  for  acceleration  response  due  to  unit  gust  velocity  normal  to  each 
aerodynamic  reference  panel,  with  the  gust  phase  taken  at  the  origin.  It  is 
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therefore  recognized  as  the  matrix  of  frequency  response  functions  for 
acceleration, 


[H]  =  -u)2(q/V)  [RJ  [B] -1  [G]  (3.5) 

(ACCELERATION  RESPONSE  FUNCTIONS) 

3.  3  ELASTIC  LOADS 

_ 2 

Solving  for  the  elastic  loads  [M]  C'to  >]  {£}  in  equation  (3.2)  and  substitut¬ 
ing  equation  (3. 3)  for  {£}  yields 

(q/V)  [M]  fall  [B]_1[G]  {W}  =  (q/V)  ( [G]  +  [E]  [B]  _1  [G]  )  hTiQx>J  {W} 

(3.6) 

where  [E]  =  [M]  (uj2  -  }  +  q  [A]  +  q[C][F][S] 


The  matrix  coefficients  of  (W)  on  each  side  of  equation  (3.6)  contain  the 
influence  coefficients  for  elastic  load  response  due  to  unit  gust  velocity 
normal  to  each  aerodynamic  reference  panel,  with  the  gust  phase  taken 
at  the  origin.  They  represent  alternative  forms  of  the  frequency  response 
function  for  elastic  loads  and  correspond  to  the  mode  displacement  and 
mode  acceleration  methods,  respectively  (Reference  9).  Accordingly, 

[H]  =  (q/V)  [M]  [B] _1  [G]  (3.7) 

(MODE  DISPLACEMENT  METHOD) 

[H]  =  (q/V)  (  [G]  +  [EJ  [B]  ”1[G) )  (3.8) 

(MODE  ACCELERATION  METHOD) 
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The  mode  displacement  method  is  sometimes  preferred  because  it 
requires  only  the  generalized  form  of  the  aerodynamic  load  coefficients, 
whereas  the  mode  acceleration  method  requires  both  forms.  However, 
to  =  0  for  the  rigid  body  degrees  of  freedom,  so  that  according  to 
equation  (3.7),  the  mode  displacement  method  yields  no  loads  at  all  until 
flexible  modes  are  included.  Since  the  mode  acceleration  method  has  a 
head  start,  it  tends  to  be  closer  to  convergence  for  a  given  number  of 
included  modes.  The  two  methods  ultimately  approach  the  same  limit, 
however. 
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SECTION  IV 


THREE-DIMENSIONAL  GUST  RESPONSE 

The  gust  velocity  cross  spectra  and  frequency  response  functions 
developed  in  Sections  II  and  III,  respectively,  may  be  combined  to  yield 
the  various  dynamic  response  quantities  required  for  gust  responsi 
analysis.  These  include  the  response  power  spectra,  from  which  are 
derived  the  transfer  functions,  RMS  responses  and  zero  crossing  rates. 
Response  cross  spectra  with  respect  to  the  three  probe -measured 
cartesian  gust  components  may  also  be  computed  from  the  analytical 
model  if  corresponding  cross  transfer  functions  and  coherence  functions 
are  required.  These  are  us  -ul  for  verifying  the  analytical  model  when 
equivalent  dynamic  response  flight  test  data  is  available  for  comparison. 

The  aircraft  is  assumed  to  have  lateral  symmetry  to  allow  the  symmetric 

> 

and  antisymmetric  responses  to  be  computed  separately  and  superimposed. 

A  discussion  of  3-d  vertical  and  lateral  gust  response  procedure  is  speci¬ 
fically  included,  since  it  corresponds  to  the  case  treated  by  the  computer 
program  described  in  Paragraph  5.  3.  This  case  is  also  shown  to  reduce  to 
the  conventional  1-d  gust  response  formulation  in  the  small  aircraft  limit 
of  Paragraph  2.  5. 

4.  1  LATERAL  SYMMETRY 

If  the  aircraft  has  lateral  symmetry,  then  the  symmetric  and  antisym¬ 
metric  responses  may  be  computed  separately  and  superimposed.  Quanti¬ 
ties  pertaining  to  the  symmetric  and  antieymmetric  response  calculations 
will  be  distinguished  by  affixing  +  and  -  superscripts,  respectively. 

In  particular,  each  row  of  the  matrix  [H]  in  Section  III  corresponds 
to  a  given  response  and  contains  the  elements  Hm(f2)  ,  where  m  =  1, .  .  .  ,2N  , 
and  N  is  the  number  of  aerodynamic  reference  panels  on  one  side  of  the 
aircraft.  In  the  symmetry-restricted  case,  however,  the  matrices  of  the 
symmetric  and  antisymmetric  frequency  response  functions  are  written 
IH]  ±  ,  and  a  row  of  elements  is  given  by 
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(4.1) 


H±  (ft)  =  H  (Q)  ±  H  ,(fl)  (m  =  1,...,N;  in'  =  m  +  N) 
m  m  m 

Here,  each  column  corresponds  to  panel  pair  m  ,  consisting  of  aerodynamic 
reference  panel  m  and  its  image  m'  =  m  +  N  reflected  in  the  plane  of 
symmetry. 

* 

The  single  panel  frequency  response  functions  may  be  recovered  by 
alternately  adding  and  subtracting  the  two  equations  (4.  1)  to  yield 


=  K(n)  + 


(4.2) 


Since  the  gust  phase  of  each  frequency  response  function  is  taken  at  the 
origin,  the  longitudinal  separation  may  be  removed  from  the  correspond¬ 
ing  gust  velocity  cross  spectra.  Accordingly,  let 


(ft)  =  *(n_'n«'Enm'ft)  (m  =  1 , . . .  ,  2N;  n  =  1,...,2N)  (4.3) 

mn  — m  — n  Man 

This  expression  defines  the  cross  spectrum  between  gust  velocity  components 
normal  to  panels  m  and  n  .  Longitudinal  separation  between  the  panels  is 
ignored  by  the  use  of  +  Zmn— 3  •  The  rem°val  of  the  longi¬ 

tudinal  separation  has  no  effect  upon  the  original  cross  spectrum  except  to 

-ifix 

eliminate  the  phase  factor  e  11111 ,  where  x _ =  x  -  x  . 

r  mn  n  m 
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Symmetry-restricted  cross  spectra  between  panel  pairs  corresponding 
to  the  frequency  response  functions  defined  by  equation  (4.  1)  are  given  by 


**  (fl)  =  («  (A)  ±  4  .  (fi)}/2 

nrn  '  mil  m  n  * 


(m  =  1 , . . .  ,  N ;  n  =  1 , .  . .  ,H ;  m '  =  m  +  N) 


(4.4) 


The  gust  velocity  cross  spectra  between  single  panels  on  the  same  side 
and  opposite  sides  of  the  aircraft,  respectively,  are  recovered  from  the  sum 
and  difference  of  equations  (4.4).  Thus, 


4 

mn 


(fi) 


4+  (Q) 

mn 


+  *  («) 

mn 


V„<a>  -  <n(a>  -  \m<a> 


Symmetry  considerations  also  imply  that 


(4.5) 


%•<“>  =  Vn<a>'  Vn'(a>  ’  V(a>'  =  V(!I)  <4‘6> 


4.  Z  RESPONSE  POWER  SPECTRA 

The  power  spectral  density  of  an  aircraft  response  to  gust  is  given  by 
the  quadratic  form  (References  1-4), 


2N  2N 


<M«>  =  l  l  H  (fi)  4mn(n)Hn(il) 
m=l  n=l 


(4.7) 


N  N 


=  y  y  (H  $  H  +  H  4  ,11  ,  +  H  ,4>  ,  H  +  H  ,  $  ,  .  H  ,  ) 

m=l  n=l  mnrnn  m  mn'  n'  m'm'nn  m'  m'n'  n,; 


Z4 


where  the  reduced  frequency  argument  is  suppressed  for  brevity,  and  the 
single  term  is  expanded  into  four  by  using  a  primed  subscript  to  denote 
the  image  panel.  Substitution  of  equations  (4.2),  (4.  5)  and  (4.6)  into  (4.7) 
yields  after  some  manipulation, 


where 


4>(fl)  =  <P  (fl)  +  <T(ft) 

+  N  N  + 

4>  (fi)  =  I  l  H~  (0)*"  (fi)H‘(fl) 

m  mn  n 

m-1  n=l 

r>  (  \  ±  •  2  i  ,  +*  +  +  ^  +  +  . 

“.'■l  m1  nun  £.  m  lhn  n  n  nn  * 
m=l  n=l 


(4.8) 


*  J,  {  |Hm|2tL+  2t  (»e<»rHn,Re(C>  -  Im(Hm*Hn,Il"(t™,»)} 


Notice  that  the  power  spectra  <JT  (fi)  of  the  symmetric  and  antisymmetric 
responses  are  computed  separately  and  summed  to  obtain  the  total  power 
spectrum  <j>  (Q)  . 

The  computational  effort  required  to  generate  the  N  (N  +  1)  gust 
velocity  cross  spectra  required  in  the  two  response  expressions  is  usually 
negligible  compared  to  the  rest  of  the  calculation.  Furthermore,  gust 
input  panels  having  the  same  transverse  coordinates  and  spatial  orientation 
may  be  combined  by  summing  the  frequency  response  functions  which  are 
referenced  to  them.  This  can  result  in  a  significant  reduction  in  the  effec¬ 
tive  number  of  panels.  For  example,  a  fuselage  whose  constituent  panels 
are  lined  up  in  a  longitudinal  strip  may  be  treated  as  a  single  panel. 

4.3  TRANSFER  FUNCTIONS 

Since  symmetric  and  antisymmetric  responses  are  predominantly  caused 
by  vertical  and  lateral  gust,  respectively,  it  is  appropriate  to  define  corres¬ 
ponding  transfer  functions 
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F 


t 


T+(il)  =  4+  (fl)/*3(f2) 
T_(fl)  =  4“(fl)/*2(ft) 


(4.9) 


This  allows  equation  (4.8)  tc  be  cast  into  the  alternative  form 


<MQ)  =  t  2+(fi)$3(fi)  +  T2"(n)$2(n) 


(4.  10) 


The  response  transfer  functions  are  usually  quite  insensitive  to  the  turbu¬ 
lence  field  parameters,  particularly  if  the  scale  of  turbulence  is  much 
larger  than  the  transverse  dimensions  of  the  aircraft.  This  suggests  that 
the  analytical  transfer  functions  may  be  used  with  measured  gust  power 
spectra  in  equation  (4.  10).  Thus,  if  the  aircraft  is  equipped  to  measure 
both  gust  component  and  response  time  histories,  response  power  spectra 
may  be  extracted  from  dynamic  response  flight  test  data  for  direct  com¬ 
parison  with  corresponding  analytical  power  spectra,  thereby  providing 
a  means  of  substantiating  the  analytical  model.  This  procedure  is 
followed  in  the  comparison  study  described  in  Section  V. 

4.  4  RESPONSE  PARAMETERS 


The  RMS  load  amplitude  o  ,  and  the  characteristic  frequency  (or  zero¬ 
crossing  rate)  Nq  represent  the  significant  dynamic  response  parameters 
applicable  to  fatigue  design  and  limit  load  calculations.  These  quantities 
may  be  obtained  directly  from  the  corresponding  response  power  spectra 
by  the  following  expressions. 


o  = 


(Mn)dn 


>) 


1/2 


No  "  2ia  (f  *(Q)fi2dn 


1/2 


(4.11) 
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t. 


When  calculating  RMS  load  amplitudes,  it  is  customary  to  employ  a  finite 
upper  limit  on  the  integral,  since  the  power  spectrum  of  a  typical  gust 
load  tends  to  diminish  rapidly  with  frequency.  However,  the  integral  used 
to  calculate  the  zero  crossing  rate  is  usually  quite  sensitive  to  cutoff 
frequency  because  of  the  squared  frequency  factor  present  in  the  integrand. 
As  a  result,  the  high  frequency  tail  of  the  response  power  spectrum  may 
significantly  affect  the  zero  crossing  rate,  even  though  its  contribution  to 
the  RMS  amplitude  is  negligible.  To  eliminate  this  high  frequency,  low 
amplitude  influence,  the  upper  limit  of  the  zero  crossing  rate  integral  is 
typically  set  to  a  frequency  corresponding  to  0.99 o  (Reference  10). 

4.  5  RESPONSE  CROSS  SPECTRA 

If  a  gust  probe  is  mounted  on  the  aircraft  at  point  (xo,yo,ZQ)  *  then 
in  analogy  with  equation  (4.  3),  let 


=  * (e-j'An'EQn'ft)  <j  =  1,2,3;  n  =  1,...,2N)  (4.12) 

This  expression  defines  the  cross  spectrum  between  cartesian  gust  velocity 
component  j  measured  at  the  probe  and  the  gust  velocity  component 
normal  to  panel  n  .  Longitudinal  separation  between  the  probe  and  the 
panel  is  ignored  by  the  use  of  £Qn  =  yonf-2  +  Zon— 3  * 

If  the  probe  is  mounted  in  the  plane  cf  symmetry  (y  =  0)  ,  then  the 
first  of  equations  (4.6)  and  further  symmetry  considerations  yield 


$  . 
Dn 


I 


(n) 


Vn(fi) 


(j  =  1,3) 

-vn) 

(j  =  2) 

(4.13) 
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The  cross  spectrum  between  a  probe-measured  gust  velocity  com¬ 
ponent  and  an  aircraft  response  is  then  given  by  (References  1-4) 


i  fix  2N 

=  e  °  [  *.  (f2)Hn(n) 

J  n=l  J 


(4.  14) 


iHx„  N 

=  e  Q  l 

n-1 


($  .  H  + 
jn  n 


) 


The  above  cross  spectrum  expression  may  be  compared  with  equation  (4.7) 
for  the  response  power  spectrum.  Notice  that  the  summation  is  single,  and 
that  a  phase  factor  is  required  to  shift  the  gust  phase  of  the  probe  measure¬ 
ment  to  the  origin  of  the  aircraft  coordinate  system.  Substitution  of 
equations  (4.2)  and  (4.  13)  reduces  (4.  14)  to 


<p.  (ft) 


N 

°  I 

n=l 


4>  . 
Dn 


(ft)  x 


H*(ft) 


Hn(ft) 


(j  =  1,3) 
(j  =  2) 


(4.15) 


4.  6  CROSS  TRANSFER  FUNCTIONS 

Rather  than  to  compare  corresponding  response  cross  spectra  obtained 
from  dynamic  response  test  and  analysis,  it  is  somewhat  more  meaningful 
to  compare  the  corresponding  cross  transfer  functions,  given  by 


0^  (n)  =  (p.  (n) 


(4.  16) 


Physically,  this  function  represents  the  complex  amplitude  ratio  between 
the  statistically  correlated  portion  of  the  response  and  the  corresponding 
gust  velocity  component  j  measured  at  the  probe. 
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Another  useful  function  for  comparison  purposes  is  the  fraction  of 
the  mean  square  response  amplitude  that  is  statistically  correlated  with 
gust  velocity  component  j  measured  at  the  probe.  This  quantity  is 
known  as  the  coherence  function  and  is  written 


. 2 (a)  =  |  <f>j  («)  1 2/<j>  (n)*j  (ft) 


(4.17) 


4.  7  VERTICAL  AND  LATERAL  GUST 


The  vertical  and  lateral  gust  response  case  is  of  particular  interest 
because  response  to  longitudinal  gust  is  usually  omitted  from  the  analytical 
representation  of  the  aircraft.  This  case  will  be  considered  in  detail, 
since  it  is  the  formulation  used  in  the  computer  program  and  flight  test 
comparison  described  in  Section  V,  and  it  provides  a  convenient  recapitula¬ 
tion  of  the  preceding  development. 

If  longitudinal  gust  is  ignored,  then  by  the  procedure  described  in 
Paragraph  2.  7,  the  gust  velocity  cross  spectrum  given  by  equation  (2.  1) 
reduces  to  a  real  function  except  for  the  phase  factor  for  gust  penetration. 
Consequently,  if  this  result  is  substituted  into  equation  (4.  8)j  the  term 
which  contains  the  imaginary  part  of  the  gust  velocity  cross  spectrum 
vanishes.  If  this  expression  is  in  turn  substituted  into  equations  (4.9)  and 
the  definition  of  the  gust  coherence  given  by  equation  (2.4)  is  applied,  then 
the  squared  transfer  functions  become 


T2±(Q)  =  l  ( |  |  2  $1+2  l  Re(H*V)Re($*)) 

v '  m'  mm  m  n  mn  J 

m=i  n— 1 


(4. 18) 


where 


$*(«)  =  f$  (fl)  ±  $  ,  (fl) 1/2 

mn  ^  mn  m'n 
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v  =  v  -  v  «  z  =z-z,p  =/v  +z 
ymn  mn  n  m  rmn  mn  mn 


= 


^m'  ~  ~^m  '  zm'  “  Zm  '  nm'2  nm2  '  nm'3  nm3 


n  -=  -sin  y  ,  n„,  =  cos  y 

m2  'm  m3  'm 


(3-D  VERTICAL  AND  LATERAL  GUST) 


Notice  that  the  vertical  and  lateral  components  of  the  unit  normal  vector 
are  expressed  in  terms  of  the  dihedral  angle  of  the  reference  panel. 

The  cross  transfer  functions  of  equation  (4.  16)  become  by  a  similar 
procedure, 


C,(ft)  =  e 


iftx  N 

°  l 

n=l 


3n 


H^(n) 

(j  =  3) 

n 

Hn(n) 

(j  =  2) 

where 


=  VPon'“)n„j  +  ®N(P0n'n)(Vn2  +  zonnn3> 
(3-D  VERTICAL  AND  LATERAL  GUST) 


(4.  19) 


fzon^Pon  ^ 


^n^on 


(j  = 
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Reference  to  Paragraph  2.4  shows  th<it  the  planar  and  nonplanar 
coherences  in  equations  (4.  18)  and  (4.  19)  may  be  written  in  terms  of 
their  large  scale  limit  forms. 


where 


and 


$p(p,n>  =  *p(p)  -  *N(ii)/(3/8  +  <2) 

(4.  20) 

$N(p,Q)  =  *N(u)(l  +  k2) / (3/8  +  K2) 

*p(y)  =  0 . 9944  y5/6  K5/g (y) 
iN(y)  =  0.9944  (3/8) y13 /6  K1/g(y) 

<  =  1.339LA,  y  =  (p/1 . 339L)  +  k1 


Anisotropic  turbulence  as  discussed  in  Paragraph  2.8  is  introduced 
by  setting 

L  =  L3  (4.21) 


by  /$2  (J2)  /$3  (£2)  in  the  symmetric  case 
by  /¥^mvr2  m  in  the  antisymmetric  case 

(ANISOTROPIC  TURBULENCE) 


and  by  multiplying: 


n  -  and  n  „ 
m2  n2 


n  and  n  _ 
m3  n3 
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If  1*2  and  Lj  are  arbitrarily  set  to  very  large  values,  the  result 
will  correspond  to  the  large  scale  limit  of  Paragraph  2.  3  in  which  the 
scale  of  turbulence  disappears  as  a  parameter,  and  equations  (4,20) 
are  replaced  by 


*p(p#ft>  =  4>p(p*<> 

(LARGE  SCALE  LIMIT)  (4.22) 

$N(p,fi)  =  $N(pfi) 


If  the  lateral  and  vertical  coordinates  of  the  gust  probe  and  the  aero¬ 
dynamic  reference  panels  are  arbitrarily  set  to  zero,  the  result  will 
correspond  to  the  small  aircraft  limit  of  Paragraph  2.4  in  which  the  scale 
of  turbulence  again  disappears  as  a  parameter,  and  equations  (4.20)  are 
reduced  to 


$p(p,f2)  =  1 

(SMALL  AIRCRAFT  LIMIT)  (4.23) 

*N(p,n>  =  0 


This  is  identical  to  the  1-d  turbulence  case  discussed  in  Paragraph  2.6. 
Substituting  this  result  into  equations  (4.  18)  and  (4.  19)  yields  the  conven¬ 
tional  1-d  vertical  and  lateral  gust  formulation, 


T+(ft) 

T"(«) 


(C3  (£2) 
|C2(fi) 
C.  (ft) 


(1-D  VERTICAL  AND  LATERAL  GUST) 


iftx  N 

e  0  7  n  .  x 
n=l  n=> 


fH* (ft)  (j  =  3) 


H~ (ft)  (j  =  2) 


(4.24) 
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where 


SECTION  V 


EVALUATION  OF  THE  3-D  MODEL 

The  three-dimensional  gust  response  methods  described  in  the  pre¬ 
ceding  sections  were  incorporated  into  a  dynamic  response  computer  pro¬ 
gram  which  was  utilized  to  obtain  a  comparison  between  1-D  and  3-D 
analytical  results.  These  results  were  substantiated  by  further  comparison 
with  flight-measured  data  from  the  C-5A  aircraft.  Because  of  its  size,  con¬ 
figuration,  and  subsonic  speed,  this  aircraft  provides  a  unique  opportunity 
to  demonstrate  and  evaluate  the  effects  of  three-dimensional  turbulence. 

5.  1  AIRCRAFT  DESCRIPTION 

The  C-5A  is  a  land  based  cargo  transport  with  a  four  turbofan  jet 
engine,  high-wing  monoplane  configuration.  The  empennage  is  a  canti¬ 
lever  Tee  arrangement  of  a  vertical  stabilizer  and  a  variable  incidence 
horizontal  stabilizer. 

The  total  wing  span  is  2630  inches;  the  basic  fuselage  length  is  2767 
inches;  and  the  horizontal  stabilizer  span  is  812  inches.  The  four  engines 
are  attached  to  wing  pylons  and  are  located  two  on  a  side  at  wing  butt  lines 
476  and  743.  Sweep  at  the  local  quarter  chord  is  25.  0°  for  the  wing,  25.  5° 
for  the  horizontal  stabilizer,  and  34.  9°  for  the  vertical  stabilizer.  Equipped 
weight  empty  of  the  aircraft  is  323,  904  pounds. 

5.  2  COMPUTATIONAL  MODEL 

The  flight  test  event  used  in  the  comparison  is  characterized  by  a  set 
of  conditions  which  are  incorporated  into  the  corresponding  theoretical 
case.  These  consist  of  altitude,  gross  weight,  fuel  weight,  cargo  weight, 
Mach  number  and  airspeed. 
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The  theoretical  model  represents  the  aircraft  structure  by  the  use  of 
six  rigid  body  and  the  first  fifteen  symmetric  and  fifteen  antisymmetric 
flexible  modes  of  vibration,  based  on  a  system  of  68  lumped  masses  con¬ 
nected  by  flexible  shafts  having  the  elastic  and  geometric  properties  of  the 
intervening  structure. 

The  analytical  aerodynamic  configuration  of  the  C-5A  analytical  model 
is  shown  in  Figure  4.  It  is  constructed  from  a  series  of  plane  trapezoidal 
surfaces  which  are  divided  into  span  and  chord  segments  to  define  subpanels. 
Doublet  lattice  methods  (References  11  and  12)  are  applied  to  the  resulting 
subpanel  distribution  to  furnish  aerodynamic  influence  coefficients  at  uni¬ 
formly  spaced  frequencies,  using  the  Mach  number  associated  with  the 
flight  test  event. 

The  resulting  mass  distribution,  modal  characteristics  and  aerodynamic 
data  are  employed  to  generate  the  inertial  and  aerodynamic  load  coefficients. 
The  inertial  load  coefficients  are  taken  with  respect  to  the  lumped  mabses, 
whereas  the  aerodynamic  load  coefficients  are  computed  relative  to  specified 
aerodynamic  reference  panels,  each  composed  of  one  or  more  of  the  doublet 
lattice  subpanels.  The  model  employed  in  this  study  contains  43  subpanels, 
which  are  combined  to  furnish  29  reference  panels. 

5.  3  COMPUTER  PROGRAM 

The  computer  program  developed  for  the  theoretical  gust  response  cal¬ 
culation  accounts  for  simultaneous  three-dimensional  vertical  and  lateral 
gust  under  anisotropic  turbulence  conditions.  The  asymptotic  planar  and 
nonplanar  coherences,  stored  in  tabular  form,  are  used  to  compute  the 
required  gust  coherences  by  correcting  for  nonzero  p/L  according  to  equa¬ 
tions  (4.20).  Frequency  response  functions  for  acceleration  are  obtained 
by  equation  (3.5).  Frequency  response  functions  for  loads  maybe  calculated 
by  either  the  mode  acceleration  of  the  mode  displacement  methods,  corres¬ 
ponding  to  equations  (3.7)  and  (3.8),  respectively.  In  the  latter  case  only 
the  generalized  aerodynamic  load  coefficients  are  required.  If  a  one-dimen- 


Figure  4. 


Aerodynamic  Configuration  for  the  C-5A 
Analytical  Model 
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sional  turbulence  field  is  requested,  the  program  furnishes  it  by  setting  to 
zero  the  vertical  coordinate  of  the  probe  and  the  vertical  and  lateral  coordi¬ 
nates  of  the  reference  panels. 

The  initial  output  data  consists  of  response  transfer  functions  and  cross 
transfer  functions,  given  by  equations  (4.  18)  and  (4.  19)  respectively. 

Response  power  spectra,  response  parameters  and  coherence  functions  are 
then  calculated  according  to  equations  (4. 10  -  4. 11)  and  (4. 16  -  4.  17),  respec¬ 
tively,  using  either  measured  gust  input  spectra  or  the  vertical  and  lateral 
analytical  spectra  given  by  equations  (2.  6)  and  (2.  18). 

The  only  significant  computational  effort  required  by  the  3-d  ^ust 
response  analysis  beyond  that  of  the  conventional  1-d  calculation  is 
attributable  to  the  generation  of  the  gust  coherences.  This  amounts  to 
approximately  3  seconds  for  a  typical  run  on  the  UNIVAC  1106. 

5.4  1-D  VERSUS  3-D  RESPONSE  PARAMETERS 

Calculated  spanwise  distributions  of  RMS  loads  and  characteristic 
frequencies  for  the  wing  and  horizontal  stabilizer  are  shown  in  Figures  5 
and  6,  respectively.  The  general  reduction  in  wing  loads  which  occurs  under 
the  3-d  turbulence  assumption  is  attributable  to  spanwise  averaging  effects, 
which  are  omitted  in  the  1-d  model.  This  interpretation  is  supported  by  the 
fact  that  the  fractional  reduction  in  loads  increases  with  span,  and  the  bend¬ 
ing  moment,  which  gives  the  greatest  weight  to  the  most  outboard  load  sources, 
shows  the  greatest  overall  decline.  Such  a  result  is  expected,  since  the  out¬ 
board  load  sources  tend  to  be  most  susceptible  to  spanwise  averaging.  The 
reduction  in  wing  vertical  bending  in  the  3-d  case  ranges  from  about  7%  at 
the  wing  root  to  about  13%  further  outboard.  Reduction  in  vertical  shear 
ranges  from  3%  inboard  to  14%  outboard,  and  torsion  from  1%  to  5%>. 

The  3-d  calculation  yields  a  totally  different  result  for  horizontal 
stabilizer  loads.  Here,  there  is  an  increase  instead  of  a  decrease  in 
vertical  bending,  and  the  fractional  change  is  a  nearly  uniform  2%>  to  3%> 
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TORSION  VERTICAL  BENDING  MOMENT  VERTICAL  SHEAR 

O'  (IN-LB  X  10  )  O'  (IN-LB  X  106)  <T  (LB  X  104) 


CHARACTERISTIC  VERTICAL  BENDiNG  MOMENT 


Figure  6.  Horizontal  Stabilizer  RMS  Incremental  Loads 
and  Characteristic  Frequencies  in  Turbulence 
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all  across  the  span.  This  is  probably  attributable  to  increased  antisym¬ 
metric  modal  excitation  caused  by  3-d  antisymmetric  vertical  gust  imping¬ 
ing  upon  the  wing. 

The  RMS  incremental  vertical  and  lateral  acceleration  results  at  the 
aircraft  center  of  gravity  (fuselage  station  1333)  are  summarized  below: 

C.G.  RMS  INCREMENTAL  ACCELERATION 
(G's)  Nq  (HZ) 


1-D 

3-D 

C  ange 

1-D 

3-D 

Change 

VERTICAL 

.0314 

.0293 

-7% 

.885 

.818 

-8% 

LATERAL 

.  0112 

.0118 

5% 

1.  574 

1.615 

3% 

5.  5  FLIGHT  TEST  COMPARISON 

The  flight  test  instrumentation  and  data  reduction  procedures  used  to 
obtain  dynamic  response  gust  flight  data  for  comparison  with  the  corres¬ 
ponding  theoretical  analysis  are  described  in  Reference  10.  Digital 
spectral  analysis  procedures  are  employed  to  obtain  flight-measured  power 
spectra  and  cross  spectra  from  which  all  of  the  necessary  flight  test  data 
may  be  computed.  Both  test  and  theoretical  procedures  are  briefly 
illustrated  in  Figure  7. 

The  spectral  method  gives  results  based  upon  the  total  output  spectra, 
without  regard  to  input  source.  The  cross-spectral  method  gives  data  based 
upon  that  part  of  the  output  that  is  statistically  coherent  with  the  input,  and 
thus  eliminates  noise  in  the  output,  regardless  of  source.  In  this  case  the 
noise  is  primarily  a  combination  of  trend  removal  errors,  uncorrelated 
pilot  inputs  and  longitudinal  gust.  Since  these  are  primarily  low  frequency 
phenomena,  they  appear  in  the  form  of  spurious  low  frequency  increments 
to  the  measured  response  power  spectra.  Because  the  coherence  properties 
of  the  turbulence  field  are  accounted  for  in  the  3-d  gust  response  analysis, 
theoretical  results  which  depend  directly  upon  the  cross  spectium  tend  to 
show  a  marked  improvement  over  the  1-d  case.  These  include  cross  trans¬ 
fer  functions  and  coherence  functions. 
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igure  7.  Procedure  for  Evaluating  Three-Dimensional  Gust  Response  Analysis  Methods 


Response  power  spectra,  cross  transfer  functions  and  coherence  func¬ 
tions  obtained  from  gust  flight  measurements  are  compared  on  the  following 
pages  with  the  corresponding  theoretical  1-d  and  3-d  gust  response  analyses. 
The  plots  include  the  vertical  and  lateral  gust  velocity  power  spectra,  verti¬ 
cal  and  lateral  accelerations  at  the  center  of  gravity,  and  loads  on  the  wing 
and  horizontal  stabilizer.  The  theoretical  loads  were  computed  by  the  mode 
acceleration  method. 

Approximately  3.7  minutes  of  recorded  data  (8925  time  history  points 
at  40  points  per  second)  were  used  to  compute  the  test  spectra  at  65  equally 
spaced  frequencies  over  a  range  of  0-5  Hz.  The  flight  test  was  conducted 
with  automatic  flight  control  systems  inactive,  and  pilot  input  at  a  minimum. 
Theoretical  and  experimental  flight  parameters  are  summarized  below. 

Altitude 
Gross  Weight 
Fuel  Weight 
Cargo  Weight 
Mach  Number 
Equivalent  Airspeed 
True  Airspeed 
RMS  Vertical  Gust  Velocity 
RMS  Lateral  Gust  Velocity 
Vertical  Scale  of  Turbulence  (Fitted) 

Lateral  Scale  of  Turbulence  (Fitted) 


2,  240  ft. 

514,  000  lb 
80,  000  lb 
110,  000  lb 
.411 

260.  9  kts 
462  ft/sec 
2.  2435  ft/sec 
2.  0196  ft/sec 
700  ft 
960  ft 
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SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

Section  V  shows  that  variations  in  incremental  loads  resulting  from 
the  application  of  a  3-d  rather  than  a  1-d  turbulence  representation  to  the 
gust  response  analysis  of  a  large  aircraft  are  far  from  uniform.  Such 
changes  cannot  be  approximated  by  adjusting  the  1-d  results  by  a  single 
empirical  correction  factor  as  has  sometimes  been  suggested.  Further¬ 
more,  it  is  difficult  to  ignore  load  changes  of  the  magnitude  ihown  in 
Section  V,  particularly  in  fatigue  design.  For  the  simpie  p'jwer  law  SN 
curve,  the  percent  change  in  fatigue  damage  is  many  times  greater 
than  the  percent  change  in  the  related  incremental  load  for  typ.cal 
aircraft  structural  alloys. 

The  3-d  turbulence  model  presented  in  this  report  is  derived  from  the 
same  turbulence  theory  as  the  conventional  1-d  model.  No  additional  assump¬ 
tions  are  introduced,  and  the  only  change  which  occurs  in  the  transition  to  the 
3-d  case  is  that  the  gust  coherence  becomes  dependent  upon  transverse  sep¬ 
aration.  Moreover,  this  dependence  is  apparently  quite  insensitive  to  the 
specific  spectral  character  of  the  turbulence  field.  For  example,  gust 
coherences  derived  from  the  von  Karman  and  Dryden  spectral  models  are 
almost  indistinguishable  according  to  Reference  4. 

In  effect,  the  3-d  model  furnishes  a  transverse  coherence  description 
which  is  mathematically  consistent  with  the  longitudinal  variation  in  the 
turbulence  field  that  is  reflected  in  the  power  spectra  of  the  three  gust  veloc¬ 
ity  components.  Since  these  spectra  are  measurable  and  their  form  is  well 
established,  the  3-d  turbulence  description  does  in  no  sense  introduce  a  weak 
link  in  the  gust  response  analysis.  On  the  contrary,  inherent  defects  in  the 
structural  and  aerodynamic  representation  of  the  aircraft  become  more 
apparent  because  they  are  no  longer  obscured  by  the  1-d  turbulence  assump¬ 
tion.  This  is  particularly  evident  in  comparisons  with  flight  test  results 
where  a  gust  probe  has  been  employed. 
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It  is  therefore  recommended  that  the  effects  of  changes  and  refinements 
in  the  analytical  model  of  the  aircraft  be  explored  within  the  context  of  the 
3-d  gust  response  analysis,  using  the  flight  test  comparison  procedures 
discussed  in  Section  V.  Development  efforts  in  the  following  areas  are 
also  suggested:  (1)  introduction  of  special  analytical  methods  to  reduce 
the  computational  effort  and  volume  of  data  which  are  ordinarily  generated 
in  gust  design  calculations,  (Z)  adaptation  of  the  current  3-d  gust  response 
method  for  time  domain  analysis  to  provide  response  time  histories,  (3) 
extension  of  the  analytical  aerodynamic  model  to  permit  the  longitudinal 
gust  component  to  be  included,  and  (4)  modification  of  the  present  3-d  gust 
response  method  to  treat  turbulence  nonstationarities  that  are  rapid  enough 
to  induce  significant  transient  responses  in  the  aircraft  structure. 
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